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Abstract. A finite-temperature effective free energy of the boundary of a quantized 
thermal Hall system is derived microscopically from the bulk two-dimensional Dirac 
fermion coupled with a gravitational field. In two spatial dimensions, the thermal Hall 
conductivity of fully gapped insulators and superconductors is quantized and given 
by the bulk Chern number, in analogy to the quantized electric Hall conductivity 
in quantum Hall systems. From the perspective of effective action functionals, two 
distinct types of the field theory have been proposed to describe the quantized thermal 
Hall effect. One of these, known as the gravitational Chern-Simons action, is a kind 
of topological field theory, and the other is a phenomenological theory relevant to 
the Streda formula. In order to solve this problem, we derive microscopically an 
effective theory that accounts for the quantized thermal Hall effect. In this paper, 
the two-dimensional Dirac fermion under a static background gravitational field is 
considered in equilibrium at a finite temperature, from which an effective boundary 
free energy functional of the gravitational field is derived. This boundary theory is 
shown to explain the quantized thermal Hall conductivity and thermal Hall current in 
the bulk by assuming the Lorentz symmetry. The bulk effective theory is consistently 
determined via the boundary effective theory. 
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1. Introduction 

Topology of the energy band strncture in insulators and superconductors emerges in 
transport phenomena, and the topological number can be detected as a transport 
coefficient [H [2]. The electric Hall conductivity is quantized and given by the Chern 
number of the occupied wave functions in two-dimensional band insulators[3l H]. The 
Chern number is nonzero when an electronic system is in time-reversal symmetry 
broken topological phases. Similarly, two-dimensional insulators and superconductors 
show the quantized thermal Hall effect in time-reversal symmetry broken topological 
phases[5]. This similarity between electric and thermal responses reflects a parallel 
between topological insulators in symmetry class A and topological superconductors in 
symmetry class D, both of which are characterized, in two dimensions, by the Chern 
number of the filled energy bands[6l [71IH]. 

Topological quantum field theories are efficient descriptions and characterizations of 
topological phases. Effective actions for the external electromagnetic and gravitational 
fields are given in terms of topological terms and attributed to quantum anomalies 
[HI [10]. They can be used to discuss the classification of topological insulators 
and superconductors in arbitrary dimensions [111 [T2l [13]. Quantized electromagnetic 
responses, including the quantum Hall effect in two dimensions, can be viewed as 
responses resulting from the Chern-Simons action functional. For the case of the 
electromagnetic field, the Chern-Simons action is given as 

= ^ld^x ( 1 ) 

in the (2-|-l)-dimensional space-time, where C is the Chern number (here and henceforth 
natural units with c = h = ks = ^ are used). On the other hand, the quantized 
thermal Hall effect has been predicted to occur in a weak pairing phase of the chiral 
p-wave superconductor in Ref. [S], which is considered to be a realization of a two- 
dimensional time-reversal symmetry broken topological superconductor. In analogy with 
the electromagnetic response, it has been claimed that the quantized thermal response 
of topological superconductors is described by the gravitational Chern-Simons action [5] 

S'^[uj] = ^ j d^xe^'^^tr ^uif,duU!p^ujf,uj^ujp^ , (2) 

where is the spin connection. In recent years, space-time curvature is widely used 
to study characteristic responses of topological materials [HI [TS] [16]. The gravitational 
Chern-Simons action is associated with the gravitational anomaly in (2-1-1) dimensions 
[9i[ini[i7], and can be microscopically derived from the (2-|-l)-dimensional massive Dirac 
fermion coupled with the background gravitational field at zero temperature [TKl [THl EU] . 
However, it has been noted that a thermal Hall current in the bulk cannot be created by 
a uniform gravitational field gradient as a response derived from the gravitational Chern- 
Simons action[2T]. Therefore, while the coefficient of the gravitational Chern-Simons 
term (“the chiral central charge”) is related to the quantized thermal Hall conductance, 
the quantized thermal transport and the gravitational Chern-Simons term appear to 
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be more remotely related then the quantized electromagnetic response and the U{1) 
Chern-Simons term. 

An inhomogeneous temperature held driving thermal transport is effectively 
realized by a gravitational potential field f through the Luttinger’s phenomenological 
argument using the Tolman-Ehrenfest relation [22] 

^dT = df. (3) 


Owing to the relation jx = Je — {F/^)jc, the thermal current Jt is identified as the 
energy current Je when the charge current jc does not contribute at zero chemical 
potential /i = 0, which is true for superconductors. A thermal current induced by a 
temperature gradient is then equivalent to an energy current induced by the space-time 
metric 

From the phenomenological point of view, an analogy between the electric and the 
thermal transport holds. The Wiedemann-Franz law connecting the electric and the 
thermal Hall conductivity for the Dirac fermion [2^ [2^ [25] has been proved to be 


TT 1 


( 4 ) 


Note that in the case of topological superconductors, the coefficient is half of that in dlj) 
since a Majorana fermion, a quasi-particle in a topological superconductor, is half of a 
complex fermion. Based on the phenomenological analogy of the Hall conductivity, an 
effective free energy for the quantized thermal Hall effect 

= KHe^^diT (5) 


of a Lorentz invariant system has been proposed as [21] 

F[0, D]oc^^- J (6) 

where D is the angular velocity of the system in a rotating frame, and k.h = Cti'^T/Qe is 
the thermal Hall conductivity for insulators whose occupied energy bands have a total 
Chern number C, and v is the Fermi velocity of the Dirac fermion. This form of the 
phenomenological effective free energy (I6|) can be realized when we consider a metric 
in a rotating frame with a gravitational potential. So far, a action or a free energy 
of the gravitational field that accounts for the quantized thermal Hall effect has not 
been derived directly from fermionic models of topological insulators and topological 
superconductors. 

In this study, an effective free energy functional of the gravitational field is derived 
microscopically from the two-dimensional massive Dirac fermion coupled with the static 
gravitational field. The effective free energy is derived by the following procedure. First, 
we consider the gapless boundary fermion, which is a manifestation of the Dirac fermion 
with nontrivial bulk energy band topology. Then an effective free energy of the boundary 
theory is calculated by the field theoretical method. We also consider the bulk thermal 
Hall current and a bulk effective field theory. 
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This paper is organized as follows. As a preliminary to the main contents of 
this study, phenomenological theory of the quantized thermal Hall effect is reviewed 
in Sec. 12.11 and a microscopic model studied in this paper is introduced in Sec. 12.21 
In Sec. 13.11 the two-dimensional bulk fermion Hamiltonian is projected onto the one¬ 
dimensional boundary to give the boundary fermion Hamiltonian, and, in Sec. 13.21 the 
boundary effective free energy is calculated from the boundary fermion. In Sec. 13.31 the 
quantized thermal Hall conductivity is derived from the boundary theory. In Sec. 14.11 
the thermal Hall current of the two-dimensional gapped bulk are deduced from the 
boundary effective free energy. In Sec. 14.21 the bulk effective free energy is conjectured 
from the expression of the bulk thermal current. We give discussion about our results 
in Sec. |5]and summarize our conclusions in Sec. |6l 


2. Preliminaries 


2.1. Phenomenology 

In this subsection, we briefly review the phenomenology of the thermal Hall effect [231 [24] 
of a Lorentz invariant system in two dimensions. Quantized thermal (energetic) Hall 
current in a two-dimensional fully gapped topological insulator and superconductor 
is phenomenologically defined in analogy with the quantized part of the electric Hall 
current in the quantum Hall system. The electric Hall current is defined by 

3c = o-ffz X E, (7) 

where z is the unit vector pointing the out-of-plane direction. Combining ([^ with the 
Streda formula[26] for the quantized electric Hall conductivity, given by the derivative 
of the magnetization with respect to the chemical potential 

dM^ 


(Th = -e- 


djj. 


and a phenomenological relation eE = —dfi, we have obtained the relation 

dM 


jc = d/ax 


dfi 


( 8 ) 


(9) 


The electric current ([H]) can be regarded as a magnetization current jc = 9 x M, if we 
assume that the magnetization varies according to variation of the chemical potential. 
Similarly, by using the Streda formula for the quantized thermal Hall conductivity |24] . 
given by the derivative of the energy magnetization Me with respect to temperature 


kh — 


dMf 
' dT '■ 


( 10 ) 


the thermal Hall current 


3e = -khz X dT (11) 

can be identified with the energy version of the magnetization current [23] Je = dx Me, 
when the energy magnetization changes through variation of temperature as 

dME 


3e = dT X 


dT 


( 12 ) 
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Note that M and Me have only ^-component. This fact implies that the thermal Hall 
current in a fully gapped system can be described in terms of the magnetization type 
energy current caused by the energy magnetization. 

Lorentz invariance of the fermionic system gives rise to a symmetry of the energy- 
momentum tensor with respect to permutations of the indices. The energy momentum 
tensor has the form of 



( 

e 


ilh 

\ 

rpfLV _ 


VTl^ 





K 




/ 


where e is the energy density, vr-^ is the momentum density, and is the stress tensor. 
The symmetry of the tensor gives the equation Je = and then it leads to an 

equality between the energy magnetization Mff = {x^Je — and the angular 

momentum as 


= v‘^L^'^12, (14) 

where Mf = and Here (• • ■) indicates evaluating an operator in thermal 

equilibrium at temperature T. Substituting the equation (IT^ and (|T4ll into the dehnition 
of the thermal Hall coefficient, we obtain 

jf dL^ dL^ 

where in the last equality of fllSp . the relation ([3]) is used. The angular momentum is 
the conjugate held of the angular velocity ff, and is given by variation of the free energy 
as = —6F/6Q^. The free energy functional for the thermal Hall effect is given by 

F[(j),Q] = — j = J (16) 

Corresponding electromagnetic free energy is the electromagnetic Chern-Simons term 
with gauge hxing oc an J d^xAoB^. The free energy represents an quantized 
electric Hall response, that is, an electric current is induced by a magnetization given 
by 


2.2. Two-dimensional Dirac fermion under gravitational field 


The Dirac fermion with the Fermi velocity v is invariant under the Lorentz 
transformation which preserves the line element 

ds^ = (ndf)^ — dx^ — dy'^. (17) 


Introducing the three coordinates x^ = {vt,x,y), a metric representing a system under 
a gravitational potential f rotating with an angular velocity D is given bv|27j 


9iJ,u 


^ l-\-2(j) Qy/v —VLx/v ^ 

VLy/v —1 0 

—Vtx/v 0 —1 y 


(18) 
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up to linear order in f and The effective field theory addressed in this paper is the 
free energy functional of the metric that results from the two-dimensional Dirac fermion 
coupled with the gravitational field. In the following, we proceed with a general metric 
form, and the metric flT 8 |) is substituted into the expression when the Tolman-Ehrenfest 
relation (I3]) is used to mimic temperature gradient. 

The two-dimensional Dirac fermion Hamiltonian coupled with the gravitational 
field is separated into the flat and the remaining parts as follows. Deviation of a curved 
space-time metric from the flat Minkowski space-time is defined by 

T h^i/. (19) 

We use the sign convention rj^y = diag(-f 1, —1, —1). The metrics and are related 
by a triad field e"^ with the equation 

( 20 ) 

When a deviation is small enough {h^^y <C 1), we can write the triad field in terms of 
h as ~ 5“^ -|- h"^/ 2 . Subscripts and superscripts of h^y are raised and lowered by 
7]^y and its inverse , like h°‘^ = hy^ri°''^. The inverse and the determinant of the triad 
field are given, up to linear order in h, by 

e/ (21) 

^ = ^det^ = det(e"^) ~ 1 + /i/ 2 , ( 22 ) 

where h = h^^. The inverse of the triad satisfies e°‘^ef' = <5)), and 

~ hoip- The gamma matrices on the curved space-time 7 ^ satisfy the 
anticommutation relation { 7 ^, 7 ^} = 2g>^'^. 7 ^ is related to the gamma matrices on 
the flat space-time 7 " by 

f = ( 23 ) 

where 7 " satisfies the relation { 7 “, 7 ^} = rj^^- 

The action of the (2-|-l)-dimensional Dirac fermion coupled with a gravitational 
field is 

S' = y dtd'^Xy/gfjj [in 7 '^V^ — m] (24) 

Here the covariant derivative of the spinor field is given by where 

the connection is given by = -'japuj/A, using 7 «^ = [ 7 a, 7 / 3]/2 with 7 « = 
and 




rP ') 


The local Hamiltonian Fi is defined by 

S = f dtd'^x^/g'^jJ^'y^'y^ {idt — Fi) fj. 


( 25 ) 

( 26 ) 


Then the Hamiltonian of the two-dimensional Dirac fermion coupled with a gravitational 
field is written as 


H 


d^x'ljJ^^/g'^y^'J^FL^p 
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= j ^ [iu7°7°a;o — if7°7-^ Vj + 7717 °] ^lJ (27) 

~ J [Ho+ U]'ijj. (28) 

Latin indices in (|27)) run over the spatial dimensions (j = 1, 2) and Greek indices run over 
both temporal and spatial dimensions {a = 0,1, 2). Matrix element of the unperturbed 
Hamiltonian is given by 

Ho = + 7717 °. (29) 

The perturbation Hamiltonian consists of two parts U = Hi + H 2 ] one from variation 
of the triad and metric in fl27)) . and the other from the spin connection. The matrix 
elements of the perturbation terms are, respectively, 

Hi = {h/2)Ho + i7;( V72)7V5i, (30) 

H 2 = -(m/4)7V7“7^ • (31) 

where ^ fhat contains at most linear order in 

and we have used = —oj^ya- Later the second perturbation term fl^ will be 

dropped out since it does not contribute to the hnite temperature free energy. Although 
the original Lagrangian has the Lorentz invariance, symmetries of the Hamiltonian (|29|) 
and fl30l) are reduced down to the SO(2) rotational invariance in the x^-x"^ plane. Time- 
reversal symmetry and parity symmetry are broken by the mass term. 

3. Effective theory at the boundary 

In the (2-|-l)-dimensional space-time, the Chern-Simons action is induced by the parity 
breaking mass term of the Dirac fermion. The Chern-Simons action of the U(l) gauge 
held is derived as an effective action by tracing out the Dirac fermionic degrees of freedom 
of the action of the (2-1-1 )-dimensional Dirac fermion coupled with an electromagnetic 
£eld[2Hl EHl EU]. Similarly, the gravitational Chern-Simons action, the Chern-Simons 
action of the spin connection, is derived by tracing out the fermionic degrees of freedom 
of the action of the Dirac fermion coupled with the gravitational held|18[ [191 EO]- 

When the theory of induced Chern-Simons terms is extended to a hnite temperature 
system by considering the imaginary time, the induced theory gives the partition 
function and the free energy of the gauge helds. However, if the fermionic system is 
fully gapped, temperature dependent part of the free energy is exponentially suppressed 
by the factor where jd is the inverse temperature and | 7 T 7 | is the magnitude of 

the gap [311 l22]- The ehective free energy of the quantized thermal Hall ehect will 
not be given in this way, since the quantized thermal Hall conductivity is linearly 
dependent on temperature. Therefore in order to derive the hnite-temperature ehective 
free energy for the quantized thermal Hall ehect, we, at hrst, consider the gapless 
boundary fermion, the existence of which is guaranteed by the nontrivial energy band 
topology of gapped fermionic systems. The ehective free energy of the gravitational held 
for the one-dimensional boundary modes is calculated by the held theoretical method. 
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This calculation procedure seems consistent with the fact that the boundary theory is 
legitimate to study the quantized thermal Hall effect [33] . 

For the purpose of deriving a hnite temperature effective free energy, we consider 
the finite temperature path integral. Tracing out fermionic degrees of freedom of the 
density matrix gives an effective free energy functional of the gravitational field as 

exp(-/3F[h]) = j V^jJ^V^|Jexp{-(3H[^lJ\^lJ,h]) 

= JjDet^"^(ia;„), (32) 

n 

where Qifujn) = + "Ho + W) is the temperature Green’s function with the 

Matsubara frequency for fermions ujn = {n + 1/2)2 ti/ f {n G Z), and Det is taken for 
the Hilbert space on the two-dimensional space and 2x2 matrix degrees of freedom. 
Expanding the temperature Green’s function with respect to U, we obtain 

hi^ ^ = In ^ ln(l -|- QqV() 

OO -| 

= + (33) 

1=1 ‘ 

where Qofun) = l/(—ia;„ -|- Hq) is the temperature Green’s function in the flat 
space-time. Then expansion of the effective free energy with respect to that is 
-P = EEo is siveii by 

( 34 ) 

n 

= ( 35 ) 

n 

where the identity In Det = Tr In is used, and Tr is taken for the same Hilbert space as 
Det. FW is independent of the gravitational field and gives only a constant, and the 
second order perturbation term is the focus of this paper. 

3.1. Boundary fermion 

In two-dimensional space, consider a boundary at = 0 between a gapped bulk 
at <0 with mass m and that at x^ > 0 with mass —m. The boundary is 
extended to an entire x^ space. The masses in the two gapped semi-infinite regions 
are smoothly connected by introducing an x^-dependent mass term, whose sign changes 
at the boundary (sgn[m(x^ < 0)] = —sgn[m(x^ > 0)]). Since the Ghern number at both 
side of the boundary differs by unity, the boundary hosts a single chiral boundary mode 
of the massless Dirac fermion. The unperturbed Hamiltonian fl2^ is decoupled into 
the x^- and the x^-dependent parts. The wave function of the boundary mode of the 
Hamiltonian fl2^ is a product of a plane wave of the x^-coordinate and a two-components 
spinor wave function of the x^-coordinate satisfying 

(—in7°7^(9i -I- m(x^)7°) go{x^) = 0. 


(36) 
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Formally the solution of 0361) is given by 






exp 




m(x'^] 


\s), 


(37) 


where we have used an identity ( 7 ^)^ = = — 1 , and the normalization constant 

is dehned hj Af = J = ifol'fo)- Here, |x^) is an eigenstate vector 

of the position operator in the Hilbert space on the x^-coordinate, and iV^o) is the 
boundary state vector of the Hilbert space on x^-coordinate with the spinor degrees of 
freedom. The two component spinor |s) corresponding to edge bound states satishes 
iy^ls) = sgn(m)|s), where sgn(m) indicates the sign of the mass in x^ < 0 , while the 
spinor satisfying iy^js) = —sgn(m)|s) corresponds to states that cannot be normalized. 
By using the relation y^y^ = — ie^’^^'jp, the condition on |s) can be rewritten in a 

convenient form, y°y^|s) = —sgn(m)|s). 

Next, we consider the Hamiltonian of the gapless boundary states resulting from the 
bulk Hamiltonian by projecting the Hilbert space on the x^-coordinate onto that of the 
boundary mode. The projected Hamiltonian is given by dehning the projection opeartor 
P = |t/^o)(t/^o| as P^HP. In the following, the x^-dependent part H = {^jJo\P'^'HP\^jJo) = 
{'ifo\'H\'if>o) of the projected Hamiltonian is referred to as the boundary Hamiltonian. 
The boundary Hamiltonian for the unperturbed bulk Hamiltonian fl2^ is 


no = 


{'fol'Hol'ifo) 

(s| — iny°y^( 92 |s) 


= iv sgn{m)d2. 


(38) 


Also we consider the projection of the perturbation term (l30l) onto the boundary mode. 
For convenience, the width of the boundary states in the x^-direction is tuned to be 
narrow enough so that typical length scale of the gravitational held is much longer than 
the width of the boundary states. This situation allows us to use an assumption that the 
metric depends only on x^ near the boundary. Then x^ and x^ is completely decoupled 
in the boundary Hamiltonian. The boundary Hamiltonian for the second term of the 
perturbation term (l30D is 

(^o|iw(V/2)7V^j#o) = -{'ifo\7°'y'^7^rn{x^Mo)vh^\x‘^)/2 

+ {s\^^r\s)wh^\x^)d2/2. (39) 

The hrst term in fl39l) is zero since it contains an integral 


poo 


rx^ 


/ dx 

'^m(x^) exp 

2v ^sgn(m) / 

dx' m(x' ) 

J —00 


Jo 






0 ^ 1=00 

(X exp 

2 n“^sgn(m 

/ dx' m(x' ) 


= 0 



Jo 


0^1 = —00 


The second term in fl39p is nonzero for a = 0, 2 due to the property of the two-component 
spinor |s). The boundary Hamiltonian for the perturbation term is written as 
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where C(x^) = [—sgn(m)h(x^)+ sgn(m)/i 2 ^(a;^) • 
for second perturbation term fl3Tl) 

= — (insgn(m)/4) ca. 


■hff{x^)]/2. The boundary Hamiltonian 


( 1 ) 

0 20 


+ 00 


( 1 ) 

1 12 


- (i^^/4) 


,(i) 


,(i) 


01 ^2 20 


( 42 ) 


3.2. Effective boundary free energy 

Here we calculate the effective free energy functional of the gravitational held of the 
boundary mode. Tracing out the boundary fermionic degrees of freedom, the second 
order perturbation term in fl3^ is given by 

( 43 ) 




where Qo = l/{—iUn + 'Ho) is the temperature Green’s function of the boundary fermion, 
and U = Hi -\- 'H 2 - Tracing over the Hilbert space on the x^-coordinate, the effective 
free energy becomes 

dg 


/3^bdry _ 






/2x 


(44) 


where q is the momentum of the gravitational held, and H*^^^ (g) is the one-loop integral 
of the single-component boundary fermion given by 

dp f 2.p -\- g' 

1 1 


n‘'’W = |E/ 


X —-^^-V- (45) 

-lUn + Vp -lUn + V{p + q) 

The leading contribution of the integral over g in fl44)) is of order q^, where we 
expand the integral in the powers of g since the long-wave length behavior of the 
gravitational held is of interest. Summation over the Matsubara frequencies gives 
icun + wp)”^ = (ddf{vp)/d{vp), where /(e) is the Fermi-Dirac distribution function 
at temperature T. The coefficient fj45ll of the order of g° is given by 


n(")(o) = 


f f dp d/(p) 


-P- 


(46) 


2v J 2% dp 

Here we have replaced vp by p for simplicity. Obviously H*^^) (0) = 0, since d/(p)/dp is an 
even function ofp. Furthermore, n(°)(0) = —fd/dirv can be neglected in this study since 
the free energy proportional to (0) is independent of temperature when substituted 
into (I44)) . which indicates that this term is not relevant to the quantized thermal Hall 
ehect. The only coefficient that is concerned is 

n'“>(o) = ^ / ^-/(P)IP' 


( 47 ) 
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where 9 is the Heaviside step function, and is regarded as the zero-temperature 
Fermi-Dirac distribution function. At low temperature, the distribution function in 
(H7|) can be expanded with respect to the temperature by the Sommerfeld expansion 
f(p) ~ 9{—p) — (7r^T^/6)(d(5(p)/dp) as 

.2rr2AU^\ ,7rT2 


nf«(o) ^ ^ ^ 

V J 2tt 6 dp 12n 

Then the effective free energy is reduced to a simple form as 


( 48 ) 


pibdry _ 


dx2[C((r2)]2 
dx^[—sgn(m)h -|- sgn(m)h^ 


12 n 

vrT^ 

48n 


' —CX) 
poo 


h\f. 


( 49 ) 


3.3. Bulk thermal Hall conductivity 

The energy current flowing along the boundary can be read off from the boundary 
effective free energy (H9|) . Using an equivalence between the energy current and the 
momentum density resulting from the Lorentz invariance, the boundary energy current 
is = nT°2 ~ 2v{5F^^'^^/5h2o), which results in 

2“"^ = ^[-sgn(m)h + sgn(m)h\ - h^f\. (50) 

The Streda formula for the quantized thermal Hall effect, the form of which is given 
in fUOj) . implies that the thermal Hall conductivity is given by ratio of the variation of 
the energy magnetization and that of temperature. Thus we consider a two-dimensional 
system with the boundary under a spatially homogeneous gravitational potential 0, 
and see how the bulk energy magnetization changes according to the change of the 
gravitational potential homogeneously. For this purpose, only the first term of the 
right-hand side of fl50|) is considered since this term is proportional to the gravitational 
potential (p. The energy magnetization is dehned by 

jE = 9 X Me- (51) 

Here, we assume the energy magnetization is uniformly distributed in the bulk as a 
response of the uniform bulk gravitational potential p. This assumption is same as the 
one used in the phenomenological theory, that is, the energy magnetization depends only 
on the temperature, which is noted just above (ITOl) . Since uniform energy magnetization 
leads to vanishing bulk energy current due to fl5T|) . the energy magnetization in the bulk 
can be evaluated entirely by the boundary energy current. 

The bulk energy magnetization can be derived from the boundary energy current 
by the following argument. Consider a boundary at = 0 separating the bulk region 
with mass m {x^ < 0) and outside region with mass —m {x^ > 0). The bulk energy 
magnetization at the region with mass m and that with mass —m have opposite sign 
since the energy magnetization changes its sign by parity transformation, which is same 
as the orbital magnetization in the quantum Hall system. The boundary energy current 







Finite-temperature effective boundary theory of the quantized thermal Hall effect 12 


is defined as the energy current flowing in x^-direction in the boundary region where 
the energy magnetization changes from bulk value to one in the outer region. Then the 
energy magnetization in the bulk x < 0 is given in terms of the boundary energy 
current as 


•bdry 

Je 


’ —OO 
poo 




— [Mf{x^ = +oo) — Mf{x^ = —oo)) 

2M|, 


(52) 


which results in 


ttT^ 

Me = -sgn(m)—ft. 


(53) 


Substituting h = h^^ = 20 and using the Tolman-Ehrenfest relation (|3]), the thermal 
Hall conductivity is then given by the Streda formula for the quantized thermal Hall 
effect: 


1 dMp ttT 

= sgn(m)-, 

T dd) ^ ^ ^ 12 ’ 


(64) 


which is the quantized thermal Hall conductivity for the Chern number C = sgn(m)/2. 


4. Bulk physical quantities 

In this section, we elaborate to make a connection between the boundary effective theory 
and the bulk quantities of the thermal Hall effect, that is, the thermal Hall current ([5]) 
and the bulk effective free energy ([6]). The expression of the bulk thermal Hall current is 
deduced from the boundary energy current fl50l) . and we conjecture a possible form of the 
bulk effective free energy which creates a bulk quantized thermal Hall current. It should 
be noted that while the bulk thermal Hall conductivity (15T|) is a direct consequence of the 
boundary effective free energy (IT9|) . results in this section may be regarded as somehow 
indirect consequence from the boundary theory, since it is unclear whether or not, at 
the thermal equilibrium, the thermal Hall current flows in the bulk or only along the 
boundary in a fully gapped two-dimensional insulator [32l [33l l34] . 

f.l. Bulk thermal Hall current 

The expression of the bulk energy current can be obtained by considering a two- 
dimensional system with the boundary under spatially inhomogeneous gravitational 
potential. If the boundary is considered as an isolated one-dimensional system obeying 
the free energy fH9|) . the energy current fl50l) can break the energy conservation law when 
the metric varies spatially: 7^ 0. Here, note that the boundary energy current is 

not dependent on time, since the boundary free energy fl49|) and responses derived from 
it are static. The energy conservation law is retained by including incoming energy flow 
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from the gapped bulk. The bulk effective free energy is determined to precisely describe 
this compensating bulk energy current. The incoming energy current from one side of 
the semi-inhnite space of the bulk is equal to the other side, since the mass term is 
inverted by the time-reversal operation or the parity inversion, which causes inversion 
of the thermal Hall coefficient. Thus, the bulk energy current Je near the boundary 
satisfies the following equation: 

= -0) = = +0) = (55) 

In the following, we show that only the hrst term of the right-hand side of (l50D 
contributes to in fl55D by taking into account that the bulk energy current, jf 
in ([55]), should not be dependent on details of the boundary, since it results from the 
nontrivial topology of the bulk energy bands. To see this, we examine how each term in 
fHfl]) changes according to the direction of the boundary. The system is rotated by an 
angle 6 anticlockwise to define a new boundary. New coordinate variables perpendicular 
and parallel to the boundary are introduced as 


X 


X' 


x^ cos 9 x"^ sin 9 

—x^ sin 9 -\- x"^ cos 9 


x‘ 


13 


R\(e)x\ 


(56) 


where R\{9) is the rotation matrix. Derivatives in the new coordinates are given by 

d'j = R^\9)dk, (57) 

where Rj^{9) is the inverse of R\{9) {Rj^{9) = R^j{—9), and Rj^{9)R\{9) = holds). 
The Hamiltonian fl2^ and (13(11) has SO(2) rotational symmetry, that is, introducing new 
vectors and tensors by 

(58) 


and 


h'\{x) = R\{9)h\{x), (59) 

h'\{x) = Rf{9)Rr{9)h^m{x). (60) 


the Hamiltonian is a scalar quantity with respect to the rotation R\{9): h'] x'] = 

Roff^h^x] and lAfy',h']x'] = /i;a:]. Note that 7° and other scalar quantities 

contained in the Hamiltonian are unchanged during the rotation. Then the effective 
free energy of the new boundary perpendicular to x'^ is given by 

ttT^ 


^bdry' _ 


dx'"|C'(x«)] 


/2m2 


12n 


(61) 


Here, = [—sgn{m)h{x'‘^) -|-sgn(m)/i'^2(^^^) ~ a-^d thus the boundary 

energy current flowing along the x'^-coordinate is 

^.bdry' ^ ^ sgn(m)/i'\ - h'^o]. (62) 

The second and third terms of the right-hand side of (l62l) explicitly depend on 
the boundary angle 9 through (15^]) and fl6U]) . while the first term does not. The bulk 
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energy current, which is proportional to the Chern number of the bulk energy band, 
should not be dependent on the details of the boundary, such as the boundary angle. 
Therefore the first term of the right-hand side of (j62l) can be regarded as a term related 
to the bulk energy current by fl55D . The same conclusion can be given by considering 
the number of indices in fl55ll . The right-hand side of the fl55l) must be a vector quantity 
with respect to rotation since the left-hand side is a vector quantity. This condition is 
satisfied when is a scalar quantity, which results in the first term of the right-hand 
side of ([62]). The bulk energy current corresponding to the boundary energy current 
^■bdry ^ —sgn(m)(7rT^/12)h is given by 

ttT^ 

Je = -sgn(m)^e^'aih, (63) 

the form of which is independent of the direction of the boundary. The equation 0631) 
exactly represents the quantized thermal Hall effect ([5|) by substituting h = 2(j) from 
the metric ffTSj) . identifying the energy current with the thermal current, and using the 
Tolman-Ehrenfest relation (|3D. The thermal Hall coefficient is given by 


1 9^ TtT 

= “f ^ 

which is quantized by the Chern number C = sgn(m)/2. 


(64) 


4 . 2 . Effective bulk free energy 

The bulk free energy for the quantized thermal Hall effect can be deduced so that it 
realizes the energy current fl6^ as a gravitational response, and is given by 

F[h] = -sgn{m)^^ / d^x e’^^hkodih. (65) 

A similar form of the bulk gravitational free energy has also been predicted in Ref. [35] . 
The free energy fl6^ is regarded as a static and Lorentz invariant version of one in 
Ref. [35]. 

Total energy conservation can be seen by the translation symmetry of the bulk free 
energy together with the boundary free energy in the temporal direction. Under an 
infinitesimal coordinate transformation 

x'^ x^ d- ( 66 ) 

the metric varies as 5hko = 5fcCo, where we assume Cf is a function of x^ and x^ since we 
consider a static theory in this study. The translation (166]) leaves the bulk free energy 
(l65]) on a region D unchanged except the boundary term on dD given by 

ttT^ f 

SF[h] =-sgn{m)j^ J d'^x e’^^{dkCo)idih) 

= -sgn(m)^ [ dl ■ iCodh). 

JdD 


( 67 ) 
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On the other hand, by dehning the boundary parallel to and the bulk in < 0, the 
translation fl66|) makes a change in the boundary free energy as 

d = d ^-sgn(m)^^ J hh2o^ 

ttT^ f 

= sgn(m)^ J dx‘^Cod2h, (68) 

which precisely cancels the boundary term in fl67ll . Therefore the total free energy is 
invariant under the translation fl66ll . and it indicates the energy of the bulk and the 
boundary is totally conserved. The boundary free energy considered here is a 

part of that represents the scalar boundary energy current shown above (16^ . while 
the other terms are irrelevant to the variation (l66|) {h^, hh'^ 2 ^ (^^ 2 )^) ^^e bulk 

thermal Hall current ((h^Q)^ and Furthermore the boundary free energy in the 

parenthesis of the left-hand side of fl68l) is a half of it, since the other half corresponds 
to the other side of the bulk {x^ > 0). 

In the presence of a gravitational potential and a rotation (lT8ll . h and hko in (l65l) are 
given by h = 20 and h^Q = {flx'^/v, —flx^/v). Therefore the integrand of fl65|) together 
with that of half of the boundary free energy is e^^hdkhio = —dfVL/v. The effective free 
energy under the metric ffTSj) is given by 

ttT^ f 

F[0, H] + H]/2 =-sgn(m)^^ J d^x 

= J d^xf^l. (69) 

The free energy (l6^ realizes the coupling of the gravitational potential and the angular 
momentum, which is the same form as the phenomenologically derived free energy 
presented in (|6]). However, it turns out that, from the held theoretical method, the 
coefhcient in (16^ is half of that predicted in |24] . 

5. Discussion 

Using the hnite-temperature held theoretical method, we have derived microscopically 
a boundary ehective free energy functional of the gravitational held, and show that 
this free energy accounts for the quantized thermal Hall ehect of a gapped fermionic 
system in two-dimension. This result gives a novel insight into the ehective theory of the 
quantized thermal Hall ehect that the boundary ehective theory, derived microscopically 
from the gapped fermionic theory, is sufficient to account for the bulk quantized thermal 
Hall conductivity, unlike previous studies of bulk ehective theories for the quantized 
thermal Hall ehect and the case of the quantum Hall ehect. However, connecting the 
boundary ehective theory with the expression of the bulk thermal Hall current and the 
bulk ehective theory requires further speculation. After making this connection, our 
boundary ehective theory is shown to be consistent with some of previous bulk ehective 
theories [241 [ 35 ] . 
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The method we have used to study the effective theory is valid for a system in 
thermal equilibrium. Thus the thermal current and transport coefficient derived from 
our effective theory are considered as those for equilibrium responses. Obviously, chiral 
boundary energy current in a topological insulator can flow even in thermal equilibrium 
as the boundary electric current does. Thus there is no doubt that the quantized thermal 
Hall coefficient in the bulk can be derived from our boundary theory, since it depends 
only on the boundary energy current. On the other hand, it is not yet clear whether the 
bulk thermal Hall current flows in thermal equilibrium or not. In our study, however, we 
have also given the expression of the bulk thermal Hall current within the equilibrium 
thermal response. 

6. Conclusion 

In this study, we have considered the Hamiltonian of the two-dimensional massive Dirac 
fermion under a static gravitational field which contains the gravitational potential 
mimicking temperature gradient, and the angular velocity coupling to the momentum 
density. We place this model in a semi-infinite plane with boundary, and consider 
thermal equilibrium at a finite temperature. 

At first, the bulk Hamiltonian is projected onto the boundary. Then the boundary 
free energy is calculated from the chiral gapless boundary fermion by tracing out the 
fermionic degrees of freedom and perturbatively expanding the free energy with respect 
to the metric. We have shown that the gravitational response of the boundary free 
energy at the boundary determines the bulk energy magnetization, and it explains the 
quantized thermal Hall conductivity through the Streda formula. 

Furthermore, the bulk energy current is evaluated from the boundary energy current 
to ensure that the total energy is conserved in the whole system. The resulting bulk 
energy current precisely represents the thermal Hall effect. The effective free energy 
functional (l65j) for the two-dimensional massive Dirac fermion is then deduced to realize 
the bulk energy current. The bulk effective free energy reproduces the same form as 
the phenomenologically derived effective free energy ([6]) by substituting a metric of a 
gravitational potential in a rotating frame. 
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